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eas PACKET SWITCHED NETWORK CONCEPT 

A new technique for data communications that has evolved 
over 10 years is called PACKET SWITCHING. In general, commu- 
nication networks may be conveniently divided into three 
types: CIRCUIT SWITCHING, MESSAGE SWITCHING and PACKET 
SWITCHING. Both message and packet switching uses a tech- 
nique known as store and forward transmission. 

A circuit switching network provides service by setting 
up a total path of connected lines from the source to the 
@eseination of the call. This complete circuit is set up by 
a special signaling message that threads its way through the 
network, seizing channels in the path as it proceeds. After 
the path is established, a return signal informs the source 
that data transmission may proceed, and all channels in the 
path are then used simultaneously. 

The entire path remains allocated to the transmission, 
and only when the source release the circuit will all these 
channels be returned to the available pool for use in other 
paths. Circuit switching is the common method for telephone 
systems. 

In message switching, only one channel is used at a time 
for a given transmission. The message first travels from 


its source node to the next node in its path, and when the 
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entire message is received at this node, then the next step 
in its journey is selected. If this selected channel (link) 
is busy, the message waits in a queue, and finally when the 
channel becomes free, transmission begins. Thus, the message 
"hops" from node to node through the network using only one 
channel at a time, possibly queueing at busy channels, as it 
is successively stored and forwarded through the network. 
Packet switching is basically the same as message switch- 
ing except that the messages are decomposed into small equal 
pieces called packets, each of which has a minimum length. 
These packets are numbered and addressed and make their way 
through the net independently of each other. Thus, many 
packets of the same message may be in tranmission simulta- 
neously, giving one of the main advantages of packet switching. 
With packet switching systems, information is exchanged 
in the form of short packets. A packet-switched network can 
handle several different types of traffic concurrently. These 
include HIGH-THROUGHPUT traffic, for example, the transmission 
of large data files between computers, for which accuracy and 
high average data speed are the most important performance 
requirements; LOW-DELAY traffic, for example, interactive 
communication between a person at a terminal and a remote 
computer, for which accuracy and low average message delay are 
important; and REAL-TIME traffic, for example, packetized 
speech for which the performance of circuit-switched connec- 


tion must be approached by maintaining a relatively constant 
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Gata speed, but for which extreme accuracy is not important 
Syvomg to the redundancy of the information. 

The packet switched network is designed primarily for 
computer to computer communication. It has a much more 
rapid response which matches the internal behavior of compu- 
ters and handles information in much the same way as does a 
computer. At the same time it can readily match the speed 
of attached computers to that of the terminal users, by virtue 
Seeits internal storage. 

The prime purvose of store and forward packet switching 
is to enable communications resources to be used effectively 
and in such a way that they may be shared by many users opera- 
meigetn an intermittent fashion, giving each user a rapid 
response from the communication network just at the instant 
when this is required. 

If there is need for transmitting a long continuous stream 
of data, then a circuit switched connection makes good sense. 
On the other hand, if the data flow is bursty, then some form 
of resource sharing can be used to great advantage; packet 
Switching is an effective choice here. 

Since packets are stored as they pass through switching 
Meces, 1t 1s possible to conduct speed, format and code con- 
version during the switching process. Another feature of 
packet switching is its ability to adaptively select good 
Paes fOr packet journeys as a function of the network 


congestion. 





Besides providing small network delays, packet switching 
has the desirable feature of rapidly handling small messages 
in spite of the presence of long messages that may be in 
transport at the same time, this is because of the decomposi- 
tion of long messages into packets. Another useful property 
of this decomposition is that the nodal storage requirement 
1s reduced. 

In evaluating packet switching networks, the delay, 
throughput, cost and reliability are important measures. 

Theoretical studies have been directed to the queueing and 
network flow problems in general and more specifically to such 
problems as delay analysis, route assignment, topological 
Gesign and flow control, etc. 

The topic chosen in this study is centered around optimal 
Flow and minimum capacity assignment and delay analysis for 
packet-switched networks. The networks under consideration 
have a relatively complicated structure with a large number 


of source-destination node pairs. 


B. THE ROUTING PROBLEM 
A message routing procedure is merely a decision rule 
that determines the node a message will next visit in its 
path through the network. The objective of the routing 
procedures is to transport packets on a minimum delay path 
rom source to destination. 
In discussing routing policies for networks, an important 


distinction must be made between static and adaptive policies. 
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This distinction depends on the environment in which a policy 
is designed to overate. If network topology is not subject 
to changes (due to failure, modifications, growth) and traffic 
inputs are stationary, then the optimal routing solution is a 
static solution consisting of a set of fixed paths between all 
node pairs. The traffic between each source and destination 
pair may be distributed on several paths simultaneously in 
well defined proportions, where the poroportions are fixed in 
time. 

In a real network environment, however, the topology 
changes with time and user traffic requirements tend to fluc- 
tuate more or less rapidly. To minimize delay it is then 
necessary to implement an adaptive routing policy that can 
meeet and adjust to various changes. The adaptivity of a 
policy can be measured in terms of its response time. A 
reasonable procedure is to use a periodicaliy refreshed 
static routing solution. Indeed, there is a continuum of 
solutions between the two extremes, static and adaptive, 
characterized by different response times and used for 
different applications. 

Beside their use in operational networks, static routing 
policies find an important application in the network design 
Process. During this process, for a given traffic pattern a 
prediction of the throughput and delay performance of a 
given topolocy is needed. The routing policy clearly has a 
major impact on such performance. Most routing implementa- 


tions are adaptive, and unfortunately the analysis of adaptive 
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routing policies 1S an extremely difficult task. To simplify 
the problem, the traffic pattern is usually approximated with 
a stationary pattern, and the routing policy with a static 
fmoemeing policy. 
ieeestactic Routing 
a. Routing Table Representation 

A static routing policy is represented by a set 
of routing tables, one for each node, indicating how the 
Meaeric Arriving at that node must be routed on the outgoing 
links, depending on its final destination. The routing table 
for node 1is an Nx L matrix Pi(n,k), where N is the number 
of nodes in the network and Pi(n,k) is the fraction of traffic 
directed to node k which, upon arrival at node i, is routed 
meeouga neighbor node n (Fig. 1.1). By the definition, 

Ba (j,k) = 1 
Where n is the number of neighbors of node i. 

MicmaccUdl ai StitouULIenl OLreLne incoming traffic 
to outgoing links may be done randomly using as probability 
weights the values Pi(j,k). If for any node i there is only 
one permissible outgoing link from i to any destination k, then 

Pi(j,k) = l for neighbor node j. 
O otherwise 
Such a routing policy is referred to as a single path 
routing policy, since only one path is used from any node i 
to any destination node k. In general, the optimal static 


PierigesOlWcLon 1S a multipath solution, allowing for the 
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simultaneous use of several routes in order to minimize 
delays. The routing table representation may also be used 


for adaptive policies. 





Eegunreel!.1., Routing Table. 


b. Problem Formulation 
In packet switched network, messages are seg- 
mented into small packets, and each packet travels from its 
source to its destination via a set of intermediate nodes. 
While awaiting transmission to the next node in the route, 
the packet must be stored until the link to that node is 
free. Thus, at each node there are several queues, one for 


Seem CUuLpUt channel. 
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Packet flow requirements between nodes arise at 
random times; therefore, link flows, queue length, and packet 
delay are random variables. A static routing problem can be 
defined as the problem of finding the static routing policy 
which optimizes the flow and minimizes the average time delay. 

Considering the relationship between routing policy 
and link flows, the problem can be formulated for solution of 
optimal flow and minimum capacity with respect to the given 
requirement as follows. 

given: (a) topology 
(b) requirement matrix R 
Mn iM Ze : alpha 
subject to: (a) conservation of flow 
(Diet tow — Capdeteye~ alpha 
Where alpha is a nonnegative scale factor. 

The mathematical model and solution techniques are 

discussed in detail later. 
me Scapeive Routing 

The adaptive routing problem consists of defining a 
procedure that dynamically updates routing tables according 
to changes in the network. This procedure will include the 
acquisition of the status of all neighboring nodes, which a 
packet associated with the source and destination node pair 
could visit next after being processed at the current node. 
The neighboring nodes are those with direct channels linking 


to a given node. 


ey, 





The structure of the routing table and the ways the 
information contained therein are to be used vary according 
to the individual routing procedure. For example, if single 
path routing is used, only one path links each source and 
destination node pair; therefore, only a single neighboring 
node associated with each node pair is listed in the routing 
table. 

The design of routing is a process of determining the 
Structure of routing tables and specifying the procedures of 
using them. The choice of routing depends on many factors 
such as network topology, available facilities, throughput 
and delay requirements, and so on. Adaptive routing tech- 


niques are ruled out in this study. 
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IIT. OPTIMAL FLOW AND MINIMUM CAPACITY ASSIGNMENT 


The optimal flow is the best selection of paths between 
source and destination, given the traffic requirements and 
the network Poh nea Milope ihe: definition of best paths may 
few ecdecoending On the nature of the traffic. In this study, 
the best paths will be regarded as the paths of minimum 
“distance”, where distance can be interpreted in various ways, 
for instance as time delay along the links. 

To investigate average time delay of the networks, two 
different approaches are discussed using classical optimiza- 


tion procedures. 


A. NETWORK MODELING 
In a distributed packet-switched network, packets are 
transmitted over a collection of switching nodes in tandem. 
In modeling a packet switched network, the objective is 
Meually to optimize the total flow of the networks, to 
minimize the capacity of each link, or to maximize the 
utilization of the network. The problem may be formulated 
as follows. 
given topology 
requirement matrix R 
capacity of each link (arbitrary) 
minimize: alpha 


subject to: (1) conservation constraints. 


HS, 





DF (k) ~  )FE€ (k) a) 

Sut gonrm6 incoming TADUEe 
links links at node i 

from node i to node i ie@omne Glen ic 


(cape cicty Gonstraints 
ETS) allies sq ele 

(3) positive constraints 

HS (US)! ea 
Where F“L(k) is the flow on link destined for node k. We 
Meme LOOK at each constraint further in detail. 

im Conservation Constraints 
At each node 1, the total outgoing packet flows 

destined for node k equal the incoming flows plus the input 


to be transmitted to node k. 


uae - )F&(k) =  Ri(k) 
outgoing incoming input 
links links at node 1 
from node 1 to node 1 For node k 
where 

Mme cesStination node (k = 1,2,3...N). 

N = number of nodes. 

L = number of links. 

je-—sGource node (1 = k, 1 = 1,2,3...N). 

vn oetween two nodes (C= 1,2,3...L). 

Pl(k) = flow on link destined for node k. 

Ri(k) = input destined for node k at node 1. 


Zo 





Geeecavacity Constraints 
The total flow destined for node k on a certain link 
Should be less than or equal to the capacity of the link ¢. 
Since linear programming algorithms have difficulties dealing 
with inequalities, slack variables are added to the constrains 
to absorb unused resources and thus to force equalities to 
appear. 


When slack variables are introduced, the problem in 


Standard form becomes. 


Heim alpoaa X Cl + Sx = 0 


where 
@ = Capacity of the link. 
alpha = constant to optimally minimize the link 
Cdederty and total tlow On fink © 
St = slack variables. 
3. Positive Constraints 
If we consider the network to have unidirectional 
links, the bidirectional flows between two nodes are divided 
into two one-way links in opposite direction, and the flow on 
Sem tink 1S positive. 
nic Use 2 18 
Using matrix notation, the model simply may be 


expressed as a linear programming problem. 


MEnamMa Ze : alpha 
subject to: AX = R 
me ae 
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B. MODEL PROGRAMMING 

In order to use a readily available Mathematical Program- 
ming System such as MPSIII, which can solve linear programming 
problems with up to 4000 rows and theoretically with an un- 
limited number of variables, the major programming problem 
with large scale networks is how to generate the input data 
for the MPSIII program. In other words, how do we generate 
the matrix form of the model. 

im Node to Link Incidence Matrix 

For any network processing to take place in the compu- 
ter, the network structure must be represented in some machine 
understandable form. A wide spectrum of different represen- 
tation methods are available. Most of these are based on an 
incidence matrix scheme. A matrix form 1S convenient for the 
analysis of networks since the resulting matrices are in a 
format suitable for mathematical analysis. 

The network of Fig. 2.1 will be used to illustrate the 
method. The node to link incidence matrix E(i,t) describes 
the links connected to the node such that 

owen eter link £ 15 OUtgoing from node i. 
eae lea Vs Incoming to node i. 
O if link £ is not connected to node i. 
The incidence matrix for input at node 1 destined for node k 
is a modification of the node to link incidence matrix £(i,¢) 


according to the following rules: 


OLED 





Hat, ) = | 0 for columns corresponding to links 
Outgoing from node k. 
delete the row if i=k 
Wome 12) otherwise 


Pemeeie illustration network (Fig. 2.1). 





E2(i,¢) and E3(i,#) may be represented in the same manner. 
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Mme 1S a 3X5 matrix and Ek(i,?) is a 2x5 matrix in size. 





Hecime: 22s Ellustration Network. 


2. Matrix Representation 
Matrix representation 1S a4 convenient tool for 
programming. The complete matrix representation for the 
illustration network is shown in Appendix A. 
The general model matrix AX = R may be represented 


as follows: 


A os = R 
El FY (1) Raat) 
E2 Ee ( 2 } = Ri (2) 
E3 Fe (3) Ri (3) 
EK Re Ok) Real (Oe) 
| ia EZ ie < Jee Ik - cl F¢ (k) 0 L 
SL | 
| alpha| 0 
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where 
Ek —Iii mn winecia@ent Matrix. 


Ik L x L identity matrix with zero elements 


COurpeSponaind to Zeno columns of the EK matrix. 
IT = L x L identity matrix. 


Ri (k) requirement matrix which is the inputs to be 


sent from source node i to destination node k. 


The size of the model in matrix representation is; 


A = ((N-1) x N+ L) x ((N+1) x L + 1) 
ae (Nt LE) x b+) xd 
R= ((N-1L) MeN Se Ep) Soe cl 


For example, the network consisting of 20 nodes and 
40 links has the matrix which is A = 420 x 841, X = 841 x1 
and R= 420 x 1 in size. 

3. Programs 

To generate the MPSTII input data and to solve the 
model, requires two programs which are listed at the end of 
this study (Appendix C and D). 

The first one is the data (model) generation program 
(Appendix C) written in the Fortran language. The output of 
this program uses special notation to designate the flow of 
the link connecting node pairs, destination node of the flow, 
flow variables, slack variables and input at the node. These 
notations are composed of a letter which is one of C, L, X or 


Seanad ¢ numerical digits. 
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iim oearietmec mii, X., ana Sc imply conservation 
Semstralnts, Capacity constraints, flow variables, and slack 
variables respectively. The rest of the six numbers following 
@aeeet, XxX and Sf indicates link connecting node pair with 
first 4 digits and destination node or slack variable number 
Wieser Last 2 digits. FOE example, the notation C1000201 indi- 
cates the conservation constraint row which is the input at 
node 2 destined for node 1, L1020100 the link capacity row for 
the link connecting between node 2 and node 1, X1020101 the 
flow variable on the link connecting between node 2 and node 
1 destined for node 1 and S1000001 the first slack variable. 

The second program is the MPSIII program developed by 
Management Science Systems, Inc. It provides elaborate con- 
trol language for the formulation of solution strategies for 
mathematical programming problems. 

Individual instructions in this language bring in 
quite elaborate sections of code designed to execute the step 
as efficiently as possible. 

Using the special notation discussed above, there is 
a limitation on the size of the problem which can be solved, 
because data names must consist of 8 characters maximum. 

This means that the maximum network size must be less than 


100 nodes, if the computer capacity allows. 


See SOLUTION TECHNIQUES 
As mentioned before, the classicial optimization proce- 


dures to solve the investigation model provide solutions 
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foci are the optimal paths and minimum link capacities for 
the network. Different optimal solutions may be expected 
depending on different methods of optimization. The two 
different approaches, each involving iteration of the linear 
programming solution procedure are discussed below, and illus- 
trated in terms of one example network problem. These 
approaches are based on the characteristics of linear pro- 
gramming problems. 

i ocuccessive Saturation Approach 

The standard procedure of the MPSIII program listed 
at the end of this study (Appendix D and E) produces an 
"optimal" solution for alpha from the output data of the data 
generation program in the first run. But the corresponding 
Flow solution is not a good solution in terms of the average 
time delay in the packet-switched network, when it is ana- 
lyzed by the Kleinrock's delay analysis model (Ref. 4). The 
problem is that only the flow thru the saturated links is 
optimum. 

In order to reduce the time delay, Iteration runs are 
attempted until all links are saturated. The detailed itera- 
tion method for the example network (Fig. 3.1) which has 5 
nodes and 6 links, each with original link capacities 10 (to 


Simplify the problem) is discussed below. 


Zz, 








Figure 3.1. Example Network. 


Dremouplt OL MPSLIt program 1s divided into row 
section and column section (Ref. 9 and Appendix D and E£). 
The objective function value (which is the same as the alpna 
value) and the links saturated with the alpha value in the 
Capacity constraint rows indicated with the first letter 
notation L can be read, and the flows on each link and the 
alpha value in column section. In the example network, for 
the first iteration the saturation level (alpha value) is 
Wez>, and the saturated links are L1010300 (link 3) and 
L1040300 (link 7). Their activity and slack activity levels 
are all zero, but their dual activity level is not (pp 48). 

For the next iteration, the changed input datas are 
that the original capacity multiplied by the saturation level 


of first run are moved from the alpha column to the right 


hand side (RHS) column. In the example, the new link capaci- 
mites 225 °(0.25 x 10) are moved to RHS (pp. 50 ), and the 
problem is solved again. Repeating this iteration procedures 


until all links are saturated and moved to RHS, the objective 
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function value becomes zero. When the objective function 
value is zero, the final flow solution is the desired solution 
of this successive saturation approach. In the example case, 
the six iterations are required to solve the problem. These 
whole procedures except programs and input datas from 2nd run 
for the example network are listed at the end of this Sia 
(Appendix D). 
2. Max~-Slack Approach 

To simplify the solution of the above successive satu- 
ration approach, another approach called the max-slack approach 
1s considered. It requires only two iteration of the program, 
but is otherwise quite Similar. Like in the successive satu- 
ration approach, after the first run, all saturated link 
Capacities are multiplied by the alpha value of the first run 
are moved to the RHS, and in addition, all slack variable 
columns are summed as the objective function in the data card 
for the second run. 

In the example network, eight slack column variables 
are added on the objective rows, and all new capacities 2.5 
(0.25 x 10) are moved to RHS. This is clearly shown on the 
Output of control language "picture", and the program, input 


data and the output are also listed at the end (Appendix E£). 


D. DELAY ANALYSIS 
The average time delay of a network is the average time 
a packet spends in the network traveling from its source to 


its destination nodes. To obtain a tractable expression for 
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average time delay T in terms of the capacity of link ¢, it 

is necessary to make simplifying assumptions. Packet lengths 
are assumed to be exponentially distributed, and are re-chosen 
With statistical independence at each node to form a Poisson 
process. 

These assumptions known as the ‘Independence assumption', 
introduced by D. Kleinrock (Ref. 4), make each queuing prob- 
lem independent. In the analysis, each queue in a packet 
network is assumed to be an M/M/1 system with arrival rate 


of Gl and to be independent of each other. 


The average delay for link 2 is given by the waiting time in 
the queue as 


ee 
uC p—Gp 


Tp = 

Where Gl is the packet traffic in the link and 1/u the 
average packet length in bits. To form a suitable average 
over all the gueuing processes, the T are weighed by Gt/r 
where r is the total input packet rate to the network. In 
this way, the total delay suffered by all packets per second 
of network operation ) Gi-T is divided by the total number 
of packets carried by me network per second. 


The total average time delay per packet through the net- 


work is then given by 


y ve ot 
PZ1 uc p - Gp BS ce 


rj 
it 
ay ae 


30 





Where Gf/u = FZ, and L is the total number of links in 

the system and 1/u is equal to the average flow on link ¢. 
It should be noted that the value of Ct (the minimum required 
Capacity) must be greater than the total of all the flows on 
link £. In other words, the value of Cl must be greater than 
the value which is the first run saturation level multiplied 
Byeeeme arbitrary capacity C for the first run. 

The value of T includes both of the waiting time and the 
service time. The waiting time is subject to the interference 
of all other traffic within the network that consists of the 
data traffic as well as the control traffic. This average 
time delay can be used to analyze the relative goodness of the 
routing strategies considered in this study. 

In deriving this basic formula, a number of factors (such 
as processing time and propagation delay) are neglected. In 
any realistic network, these variables as well as others must 
be considered in the analysis of the networks. According to 
the above formula, the average time delay depends only on the 
aggregate flow of the links. The reduction of the flow Ft of 
each link results in the reduction of the average time delay 
for the network which has constant input packet rate and link 
capacities. 

For the example network (Fig. 3.1), the arbitrary capa- 
city is 10 and the value of Ct must be greater than 2.5 
(0.25 x 10). The average time delay for the capacity FL=10, 


and the input (requirement matrix) of 2 units from node l to 


ee 





meees and 5 units from node 2 to node 3 are Tss = 2.12/r 
for the successive saturation case and Tms = 2.083/r for the 
max-slack case. To make these calculations easy, the data 
summarizing the final output run are listed in Table I. We 
see that there is not much advantage of one procedure over 
the other in terms of performance, whereas the max-slack 
solution requires only two iterations, rather than 6 for the 


successive saturation algorithm. 
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Tif. EXPERIMENTAL RESULTS AND CONCLUSIONS 


A. EAPERIMENTS AND COMPARISON OF THE RESULTS 

The experiment was conducted for two experimental networks 
to analyze both approaches, in addition to the example network. 
These networks are specified in Appendix B. The networks were 
used in the experiment with requirement matrices of 15 units 
from node l destined for node 5, 9 units from node 7 destined 
For node 2 for the 9 nodes/36 links network (Fig. B.1): and 
15 units from node 2 destined for node 12, 4 units from node 
5 destined for node 6, and 10 units from node 10 destined for 
node 1 for the 13 nodes/60 links network (Fig. B.2). Both the 
successive saturation and max-slack approaches were used for 
the successive saturation approach, the 21 iterations are 
required in this experiment. 

The results and their comparison is listed in Tables II, 
III, and IV. These are obtained by writing each flow variable 
value from the final output run on each link of the networks, 
and summing these flow variable values to obtain the total 
flow on the link (Appendix F and G). 

As shown in Tables II and III, the number of iterations 
required and the total link usage (utilization) are much 
greater for the successive saturation approach, but fewer 
large link capacities are required than in the max-slack 


approach. The average time delay depends on the capacity 
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value Cf (Table IV). When the capacity value Cf is close 
to the minimum required capacity, the time delay is less for 
the successive saturation approach than for the max-slack 
approach in both networks, and vice versa. It is significant 
for the large scale network, and for large enough link capa- 
Cities Cl, that the ene delay difference becomes negligible. 
Also, in the process of this study, it was noted that both 
approaches yield the same aggregate flows in the case of very 
Simple networks (nodes less than or equal to 4), and for com- 
pletely symmetric networks and requirement matrices (input and 
Capacity matrices). For example, the time delays are the same 
for both approaches with 2 units input from the outmost nodes 
memene directly opposite nodes respectively (1 to 13, 2 to l2, 
mmeo 9, 12 to 2 and 13 to 1) for the 13 node network, because 
all links total flow are 3/4. The detail experimental results 


are listed in Appendix H. 


fee CONCLUSIONS 

This is an elementary study of the routing design problem 
for a packet network: given a traffic requirement matrix, 
minimize the average time delay per packet, subject to finding 
a feasible flow for a network with fixed topology and link 
capacities. 

The general performance characteristics and advantages 
of the two approaches are investigated for small and simple 
Networks. This illustrates one way to solve the design 


Droblem under limited conditions. 
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TABLE IV 


TIME DELAY FOR THE CAPACITIES 
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APPENDIX A 
MATRIX REPRESENTATION FOR ILLUSTRATION NETWORK 
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Matrix Revresentation for Illustration Network. 
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APPENDIX B 


EXPERIMENTAL NETWORKS 





mucquce &.l. The ? Node/36 Link Network. 





Figure B.2. The 13 Node/60 Link Network. 
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